Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 18 Apiil 2008 (MN style file v2.2) 



Towards the optimal window for the 2MASS dipole 

Michal J. Chodorowski/* Jean-Baptiste Coiffard,^ Maciej Bilicki, 
Stephane Colombi^ and Pawel Cieciel§ig^ 

^ Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland 
^ Universite Paris XI, 91400 Orsay, France 

^ Institut d'Astrophysique de Paris, CNRS, 98 his Boulevard Arago, 75014 Paris, France 



00 ' 18 April 2008 
O 
O 

< 

00 



6 



(N 
> 

o 

o 
l> 
o 



X 



ABSTRACT 

A comparison of the 2MASS flux dipole to the CMB dipole can serve as a method to constrain 
a combination of the cosmological parameter and the luminosity bias of the 2MASS 
survey. For this constraint to be as tight as possible, it is necessary to maximize the correlation 
between the two dipoles. This can be achieved by optimizing the survey window through 
which the flux dipole is measured. Here we explicitly construct such a window for the 2MASS 
survey. The optimization in essence reduces to excluding from the calculation of the flux 
dipole galaxies brighter than some limiting magnitude A'min of the near-infrared Kg band. 
This exclusion mitigates nonlinear effects and shot noise from small scales, which decorrelate 
the 2MASS dipole from the CMB dipole. Under the assumption of negligible shot noise we 
find that the optimal value of Xmin is about five. Inclusion of shot noise shifts the optimal 
iCniin to larger values. We present an analytical formula for shot noise for the 2MASS flux 
dipole, to be used in follow-up work with 2MASS data. 

The misalignment angle between the two dipoles is a sensitive measure of their correla- 
tion: the higher the correlation, the smaller the expectation value of the angle. A minimum 
of the misalignment is thus a sign of the optimal gravity window. We model analytically the 
distribution function for the misalignment angle and show that the misalignment estimated 
by Mailer et al. is consistent with the assumed underlying model (though it is greater than 
the expectation value). We predict with about 90% confidence that the misalignment will de- 
crease if 2MASS galaxies brighter than i^min = 5 mag are excluded from the calculation of 
the flux dipole. This prediction has been indirectly confirmed by the results of Erdogdu et al. 
The measured misalignment constitutes thus an alternative way of finding the optimal value 
of -K^niin: the latter corresponds to a minimum of the former. 

Key words: methods: analytical - cosmology: large-scale structure of Universe - cosmol- 
ogy: cosmic microwave background - galaxies: general - galaxies: infrared - galaxies: Local 
Group 



1 INTRODUCTION 

The dipole anisotropy of the cosmic microwave background (CMB) 
is interpreted as a direct measure, via the Doppler shift, of the mo- 
tion of the Local Group (LG) relative to the CMB rest frame. The 
components of this motion of non-cosmological origin (the motion 
of the Sun in the Milky Way and the motion of the Milky Way in 
the LG) are known and can be subtracted (e.g., Courteau & van 
den Bergh 1999). When transformed to the barycenter of the LG, 
the motion is towards (I, b) = (273° ± 3°, 29° ± 3°), and of am- 
plitude Vi^G — 627 ± 22 km • s~^, as inferred from the first-year 
WMAP data (Bennett et al. 2003). 

The kinematic interpretation of the CMB dipole is strongly 



supported by its remarkable alignment with the dipole component 
of the large-scale galaxy distribution (often called the 'clustering 
dipole'), inferred from various all-sky surveys. In the gravitational 
instability scenario, this alignment is expected: peculiar velocities 
of galaxies are induced gravitationally and are thus strongly cou- 
pled to the large-scale matter distribution. Linear theory predicts 
the peculiar velocity of the LG, v, to be proportional to the LG pe- 
culiar acceleration, caused by the gravitational pull of surrounding 
matter inhomogeneities. Let us denote by 5 the mass density con- 
trast, 5 = q/ Qb — 1, where q is the mass density of matter and Qb 
is its average value. The clustering dipole. 



(1) 
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is a quantity proportional to the peculiar gravitational acceleration 
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(so we will call it interchangeably 'scaled gravity'), and can be esti- 
mated from a three-dimensional all-sky galaxy survey. In the linear 
regime, the relation between the velocity and the scaled gravity is 



V = //o/(an)g. 



(2) 



Here, Hg is the Hubble constant, Q^n is the cosmic matter density 
parameter and f{flrn) — ^m^' (e.g., Peebles 1980). For a spherical 
survey J d^r r/r'' — 0, hence we can write 
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(3) 



In the following we will assume that dark matter (DM) in the 
Universe is entirely locked in DM halos of luminous galaxies. 
Modelling galaxies as point particles, the observed density field is 
si^) ~ S ^i^oi'T — Ti), where 5d is Dirac's delta; rrii and ri 
are respectively the mass and the position of the i-th galaxy. Substi- 
tuting this equation into Equation ([3} yields for the scaled gravity 



\ - rrii Vi 
47r r? 



= 81,} ,7-12; 



(4) 



thus we see that the true gravitational acceleration equals to 
47rG£i5g. We will assume further that 'light traces mass', or that 
the mass-to-light ratio for galaxies is a universal constant, T. Then 
we can write 



Qb 



Qb 



(5) 



Here, Li is the luminosity of i-th galaxy and Si is its observed flux. 
Si — Li/iTTrf. In other words, since both the gravity and the flux 
fall off as distance squared, the gravitational acceleration of the LG 
is proportional to the dipole of the light distribution (i.e., the flux 
dipole) for a constant mass-to-light ratio. The sum in Equation ([5} 
is in principle over all galaxies in the Universe, while in practice 
we have at our disposal only finite, usually flux-limited, catalogs 
of galaxies. In such catalogs, lower-mass dark matter halos will be 
underrepresented by the survey galaxies. To account for this, we 
write 



N 

Qbbh ^ 



T 



(6) 



where bi, is the resulting luminosity bias and A'^ is the total number 
of galaxies in a given survey. Combining Equation l[2j with Equa- 
tion ((6) we obtain finally (Erdogdu et al. 2006; hereafter E06) 



QhhL 



N 

■ s,u 



(7) 



In the above we have used the fact that the mass-to-light ratio 
T = Qb/ Ql, where ql is the luminosity density of the Universe. 
Equation ^ shows that in the linear theory one can predict the 
LG peculiar velocity using solely an angular (two-dimensional) 
all-sky survey, bypassing the lack of radial information, i.e. dis- 
tances. Specifically, a comparison between the CMB dipole and the 
flux dipole of a given survey can yield an estimate of the parameter 

Such a comparison was first performed by Yahil, Sandage & 
Tamman (1980) using the revised Shapley-Ames catalogue and by 
Davis & Huchra (1982) using the CfA catalogue, leading to the es- 
timates of the flux dipoles that were within 30° from the CMB 
dipole. The inclusion of redshift information, usage of progres- 
sively larger redshift surveys and theoretical improvements of the 



analyses led to smaller measured values of the misalignment. In 
particular, using the IRAS 1.2 Jy survey, Strauss et al. (1992, here- 
after S92) found that the clustering dipole points around 25° away 
from the CMB dipole. Using the further completed IRAS PSCz 
survey, Schmoldt et al. 1999 (hereafter S99) obtained the cluster- 
ing dipole within 15° of the CMB dipole. A similar analysis, based 
also on the IRAS PSCz survey, performed by Rowan-Robinson et 
al. (2000), determined the misalignment angle to be around 13°. 

Two most recent analyses of the clustering dipole employed 
the Two Micron All Sky Survey (2MASS; Skrutskie et al. 1997). In 
particular, to compute the flux dipole. Mailer et al. (2003; hereafter 
M03) used the angular 2MASS extended source catalogue, with 
a limiting magnitude of Kg = 13.57. (Approximately 740,000 
galaxies covering 90% of the sky.) E06 used the Two Micron 
All Sky Redshift Survey (2MRS): approximately 23,200 2MASS 
galaxies with measured redshifts, selected from a total sample 
of about 24,800 galaxies with (extinction-corrected) magnitudes 
smaller than Ks = 11.25. 

2MASS is the first near-infrared (JHKs passbands) all-sky 
survey. While most passbands tend to be sensitive to the instanta- 
neous star formation rate, Ks passband is most sensitive to total 
stellar mass (Bell & de long 2001; Bell et al. 2003), making this 
band a better tracer of total mass. 2MASS has an effective image 
resolution of 1" and a hundred times greater sensitivity than the far- 
infrared IRAS survey. The photometric uniformity of the 2MASS 
survey is better than 4 per cent over the entire sky including the 
celestial poles (e.g., Jarrett et al. 2003). The median depth of the 
survey is 220 Mpc (Bell et al. 2003), a distance past where the 
clustering dipole has been shown to converge 

Given all these advantages of 2MASS over other all-sky 
galaxy surveys, it is perhaps surprising that the misalignment be- 
tween the CMB dipole and the 2MASS flux dipole is not smaller 
than the corresponding one for IRAS galaxies. The value obtained 
by M03 is 16°. For the 2MRS flux dipole, E06 obtained approxi- 
mately 21°[3ln this paper we aim at answering the following ques- 
tions. First: do we understand fully the origin of this misalignment? 
Second: can one do better with 2MASS, and if so, how? 

The answer to these questions is essential for optimal esti- 
mation of the parameter /3 by comparing the CMB dipole to the 
2MASS dipole. The stronger is the correlation between the two 
dipoles, the smaller are statistical errors of such an estimate. There- 
fore, the observational window, through which the 2MASS dipole 
is measured, should be adapted to obtain the best correlation pos- 
sible. The misalignment angle is a sensitive measure of this cor- 
relation: the higher the correlation, the smaller the angle. In other 
words, a minimum of the misalignment angle is a sign of the op- 
timal 2MASS window. In this paper we will formally prove these 
statements. First, we will derive the 2MASS window. Next, we will 
optimize it under the assumption of negligible shot noise. Finally, 
we will demonstrate that a minimum of the expectation value of the 
angle corresponds to minimal variance of the resulting estimate of 



^ The inclusion of galaxy redshifts in the dipole analyses allowed the es- 
timation of the convergence depth, i.e. the distance at which most of the 
clustering dipole is generated. There is a controversy whether this conver- 
gence depth is about 50 h^^ Mpc, or rather 200 h^^ Mpc (for details see 
E06). In either case, 2MASS is deep enough to provide a reliable estimate 
of the clustering dipole. (But see Basilakos & Plionis 2006.) 

E06 computed two kinds of the clustering dipole. The second one, the 
number dipole, was even more misaligned with the CMB dipole. 
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Let us enumerate possible sources of the misalignment be- 
tween the CMB dipole and an all-sky galaxy survey flux dipole. 

• M/L 7^ const. The constant mass-to-light ratio is proba- 
bly a good assumption for (almost) all galaxies when averaged 
over many galaxies of the same luminosity. For individual galax- 
ies, however, M /L is expected to have some scatter. On the other 
hand, as stated earlier, 2MASS, unlike IRAS surveys, is mainly sen- 
sitive to total stellar mass. Consequently, the mass-to-light ratio of 
2MASS galaxies is expected to have smaller scatter than that of 
IRAS galaxies. 

• Nonlinear bias. Writing Equation ((SJ we have implicitly as- 
sumed that the total flux dipole and the magnitude-limited flux 
dipole differ in the amplitude, but not in the direction. However, 
if large-scale distribution of low-mass DM halos is different from 
the distribution of high-mass halos, then the two dipoles will not be 
collinear. 

• Nonlinear dynamics. The peculiar velocity of the LG is equal 
to the temporal integral of the LG gravitational acceleration along 
the LG trajectory. Therefore, loosely speaking, while the response 
of the LG acceleration to growing nonlinearities is 'instantaneous", 
the response of the LG velocity is time-averaged and 'retarded'. As 
a result, the acceleration of the LG is more non-linear and higher in 
amplitude (in velocity units) than the velocity of the LG (Ciecieljg 
et al. 2003). What is more relevant here, at orders higher than linear 
non-local character of gravity reveals itself and tends to misalign 
the velocity vector of the LG with the vector of its acceleration. 
However, the mean misalignment angle between the velocity and 
gravity of the LG-like regions simulated in numerical experiments 
is about 8° (Davis et al. 1991, Cieciel^g et al. 2001). 

• Observational effects: shot noise, finite volume of the survey, 
and the mask (due to the zone of avoidance, ZoA). Shot noise and 
finite volume are more an issue for 2MRS, which has a median 
depth of only 60 h'^ Mpc. Still, we devote Subsection 15.11 to a 
study of shot noise of the 2MASS dipole. M03 perform two stan- 
dard treatments of the masked area: in one of them they clone the 
sky above and below the masked region; in another they fill the 
masked region with randomly chosen galaxies such that it has the 
same surface density as the unmasked area. A recent paper by Tully 
et al. (2008) puts these methods, at least partly, in question. They 
show that there lies a void in the ZoA, which they call the Local 
Void; the LG lies on its boundary. Therefore, in a part of the mask 
there is really nothing, and filling this region with faked galaxies 
leads to a systematic error of the estimate of the LG acceleration. 
However, the role of the Local Void is a recently raised issue and 
we will study it elsewhere. 

MOB and E06 notice that the misalignment is substantially re- 
duced if they remove the brightest galaxies in the catalog. M03 re- 
move all galaxies brighter than Ks — 8 mag (375 galaxies), while 
E06 remove just five the brightest. They suggest that these galaxies 
have M/L ^ const and/or non-linearly contribute to the accel- 
eration. We will study these issues here. Specifically, the outline 
of this paper is as follows. In Section (2] we will present a formal- 
ism which will allow us to model semi-analytically the distribution 
function for the misalignment angle between the CMB dipole and 
the 2MASS flux dipole. In Section[3] we will model nonlinear ef- 
fects which appear in such an analysis. The 2MASS gravity win- 
dow will be derived and optimized in Section |4] In Section [5] we 
will account for observational errors. In Section |6] we will present 
a formal proof that, under the assumption of negligible shot noise, 
our window is indeed optimal. We will also demonstrate how to op- 
timize the window in presence of shot noise. In Section |7] we will 



show the resulting distribution function for the misalignment angle. 
A summary and conclusions will be given in Section[8] 



2 ANALYTICAL DESCRIPTION OF THE 
MISALIGNMENT 

In this Section we will model theoretically the probability distribu- 
tion function (PDF) for the misalignment angle between the CMB 
dipole and the 2MASS flux dipole. The CMB dipole estimates the 
peculiar velocity of the LG, v. The 2MASS flux dipole. Equa- 
tion estimates the gravitational acceleration - more specifically, 
the scaled gravity - of the LG, g, induced by large-scale matter in- 
homogeneities traced by 2MASS galaxies. 

Letp(g, v) denote the joint PDF for the LG scaled gravity and 
peculiar velocity. It is a standard practice to approximate it by a 
multivariate Gaussian (S92; S99). Numerical simulations (Kofman 
et al. 1994, Cieciel^g et al. 2003) show that nongaussianity of fully 
nonlinear g and v is indeed small. This is not surprising since, e.g. 
gravity is an integral of density over a large volume (Eq.lTJ, so the 
central limit theorem can at least partly be applicable. 

Using statistical isotropy of g and v, their joint PDF can be 
simplified to the form (Juszkiewicz et al. 1990; Lahav, Kaiser & 
Hoffman 1990): 



-3/2 



(27r) 



■ exp 



+ y -'^PfJ'Xy 

2(1 -p2) 



(8) 



where erg and are the r.m.s. values of a single Cartesian com- 
ponent of gravity and velocity, respectively. From isotropy, — 
is ' s) /3 and cr^ = (v ■ v) /3, where {■) denote the ensemble aver- 
aging. Next, (x, y) — (g/cTg, v/ctv), and fi — cos 9 with 9 being 
the misalignment angle between g and v. Finally, p is the cross- 
correlation coefficient of Qm with Vm, where Qm (vm) denotes an 
arbitrary Cartesian component of g (v). From isotropy. 



P = 



v) 



{g2)l/2(i;2^|l/2 

Also from isotropy. 



(9) 



(10) 



where 5mn denotes the Kronecker delta. In other words, there are 
no cross-correlations between different spatial components. 

For a given all-sky galaxy survey, the LG gravity is measured 
effectively through the window of the survey, Wg (cf. Eq.[T): 



(11) 



In contrast, the LG velocity is not estimated from a velocity survey 
(i.e., from a catalog of peculiar velocities of galaxies), but measured 
directly from the dipole anisotropy of the CMB. Still, to relate it to 
theoretical quantities, we write: 



(12) 



Here i9 = — V ■ v is the (minus) velocity divergence and we as- 
sume that the velocity field is irrotational|j Thus, similarly to g, v 



^ Kelvin's circulation theorem assures that the cosmic velocity iield is 
vorticity-free as long as there is no shell crossing. N-body simulations 
(Bertschinger & Dekel 1989, Mancinelli et al. 1994, Pichon & Bemardeau 
1999) have shown that the vorticity of velocity is small in comparison to its 
divergence even in the fully nonlinear regime. 
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can be expressed as a Coulomb (Newton) integral over its source, 
i.e. the field of the velocity divergence. Here we do not assume 
that we know the latter from observations, but we know from the- 
ory its statistical relation to the density field (see this Section and 
Section [3]l. This is sufficient for our purposes in this work. Since 
V is directly measured from the CMB dipole, the effective velocity 
window, Wv, which we have introduced in Equation il2\ . is essen- 
tially unity. (Contributions from all perturbations are included. J3 
We modify slightly this form of the window to reflect the finite size 
of the LG. Following S92 and S99, we adopt 



0, r < rLG , 

1, otherwise, 



(13) 



which has a small-scale cutoff, Tlg ~ 1 Mpc. This window is 
markedly different from those appropriate for velocity surveys: the 
latter are not spherical, have complicated shapes and finite depth 
(Sarkar, Feldman & Watkins 2007). The gravity window, Wg, of 
the 2MASS survey is derived in Section|4] 

In Fourier space, relations dill ) and J12b read: 

gk = ^<5kWg(fc), (14) 



fc2 

Vk = ^ Wv(fc), 



(15) 



where the subscript k denotes the Fourier transform. The quantity 
W is related to the window W by the following equation (S92): 



W{k) 



W{r)ji{kr)dr . 



(16) 



Here and below ji represents the spherical Bessel function of first 
kind of order /. In particular. 



Wv(fc) =jo(fcri,G). 

From equations l ll4t and il5l we have 



27r2 



and 



W^{k)P{k)dk, 



Wi{k)P^{k)dk . 



(17) 



(18) 



(19) 



Here, P{k) and P^{k) are respectively the power spectrum of the 
density and the power spectrum of the velocity divergence. Defin- 
ing 

P-a{k) 



n{k) = 



we have 



Pik) 



(v ■ v) = — ^ / W^{k)TZ{k)P{k)dk . 
2^' Jo 

Furthermore, 



(20) 



(21) 



(g • v) = 2^ / Wg(fc)Wv(fc)C(fc)Py'(fc)pi/'(fc)dfc, (22) 

where C{k) is the so-called coherence function (CF), or the cor- 
relation coefficient of the Fourier components of the gravity and 
velocity fields (S92): 



The velocity of the LG is fully nonlinear and as such cannot be approx- 
imated by low-order moments of the velocity tield. In particular, vlg is 
different from the bulk velocity of a region around it. 



C{k) 



(gk ■ Vk) 



(<5ki?k) 



(|gk|2)l/2(|v^|2)l/2 ^|5^|2)l/2^|^^|2)l/2 

Hence, we obtain finally 



P = 



^^Wg{k)W^{k)C{k)'R}'^{k)P{k)dk 







Wi{k)P{k)dk 



1/2 



l'^Wd{k)K{k)P{k)dk 



1/2 



(23) 



.(24) 



Equations l |18t , | |2U and l l24t specify all the parameters (the vari- 
ances and the correlation coefficient) that determine the joint PDF 
for g and v. Equation ((Sjl, in the absence of observational errors. 
The deviation of the correlation coefficient from unity is then due 
to different windows, through which the gravity and the velocity 
of the LG are measured, and due to nonlinear effects. The latter 
are described by two functions: the CF, and the ratio of the power 
spectra (Cieciel^g & Chodorowski 2004; hereafter C04). 

The distribution for the misalignment angle can be derived 
from the joint distribution ([8). Here we are interested in the dis- 
tribution for the misalignment angle with the observed value of the 
LG velocity as a constraint. The conditional distribution function, 
p(g|v), readily results from 



p(g|v) = (2^) 



-3/2^ 



^(1 



-3/2 



exp 



(x- 



2(1 



(25) 



(Juszkiewicz et al. 1990; Lahav et al. 1990). The distribution for 
the amplitude of the LG acceleration and the cosine of the mis- 
alignment angle is p{g, fj,\v) = 27ri7^p(g|v). The distribution for 
H is obtained by marginalizing over g. 



()i\v) = 2tt dggp{e\v), 



(26) 



and the distribution for the angle itself is p{6\'v) = 
\dii/d9\ p{jj,\-v) = sin(S) p(/ijv). This yields (luszkiewicz et 
al. 1990; Lahav et al. 1990) 

p(elv) = sm{e) exp (-9^) j-^ 

+ + q^p) exp (gV') [1 + erf (g/i)] | , (27) 



where 



py 



V2(l 



(28) 



We remind that y is the amplitude of the LG peculiar velocity in 
units of the ID velocity dispersion, y = ulg/o"v. For vi^g = 627 
km ■ s^^ (Bennett et al. 2003) and the values of the cosmological 
parameters adopted here (as described in Section^, y = 2.64. This 
might suggest that the velocity of the LG is a rare event; however, 
this is on the contrary. First, one should compare the amplitude of 
the LG velocity to the 3D velocity dispersion, (Jv,3d = \/3(7v. 
Therefore, the relevant parameter here is y' = ulg /o"v,3d ~ 1.52. 
Second, the probability that a randomly chosen region will have 
velocity greater than Ulg is equal to dy'h{y'), where h{y) = 

the Maxwellian distribution. For the lower limit 
of the integral equal to 2.64, the value of the integral is 0.07 = 7%. 
However, for 1.52, its value is 0.51 = 51%. 

The misalignment of only several degrees corresponds to a 
strong coupling between g and v. In the strong coupling limit 
1 — p<Cl, sog^l. Also, /i is then close to unity. Therefore, 
in equation ([27} we can use an asymptotic formula for the error 
function. 
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erf (s) 



1 



\/7rs 



for s > f . 



(29) 



We then obtain a small-angle approximation of the distribution for 
the misalignment angle (Lahav et al. 1990): 



with 



py 



(30) 



(31) 



Thus, in the strong coupling limit the misalignment angle, given 
the velocity constraint, is Rayleigh-distributed. The parameter 6,, 
much smaller than unity (in radians), is a characteristic measure of 
the misalignment[f| The expectation value of the angle is 



(e|v) 



(32) 



Other quantities characterizing the distribution which are of 
interest here are quantiles. In our, slightly modified notation, the 
quantile 9q of a distribution p{6) is such a number, that 



p{e) de = — . 

' 100 



For the Rayleigh distribution, 
ing quantiles are 



(33) 

0. For our purposes, interest- 



fio = \/21n(10/9)6l. , and fgo = V21nl0 6l. . (34) 



3 NONLINEAR EFFECTS 

Using numerical simulations, C04 modelled the CF (Eg. 123b and 
the ratio of the power spectra (Eq. l20l (. The simulations were 
evolved from Gaussian initial conditions. As the initial power spec- 
trum of matter fluctuations, a cold dark matter (CDM) spectrum 
was adopted (as in Eq. 7 of Efstathiou, Bond & White 1992), with 
the shape parameter P = 0.19. Both the CF and the ratio of the 
power spectra were modelled as functions of the wavevector, k, 
and the amplitude of the matter fluctuations, ag. For the CF, C04 
found the following fit: 



C{k) = [l + (aofc - aik^-"" + ai/c^)2-^l , 



(35) 



with the coefficients given by the following, power-law, scaling re- 
lations in ug : 





= 4.908 < 


0.750 
^8 






= 2.663 < 


0.734 
^8 


(36) 


a2 


= 5.889 < 


0.714 
^8 





The fit was calculated for k G [0, 1] h Mpc"^ and erg G [0.1, 1], 
with the imposed constraint C{k = 0) = 1. This constraint as- 
sures that for sufficiently large, linear scales, the relation between 
the gravity and the velocity is deterministic and linear (see Eq.|2}. 
Formula i35\ is a better fit to the CF than an earlier formula of 
Chodorowski & Ciecieljg (2002), which was less accurate for low 



^ Consequently, we could approximate 6» by 2(1 — p)/y. However, 
exact expression 13 U is not more complex, while it remains valid also for 
higher-order corrections to the distribution )30t . 



values of k. Chodorowski & Ciecieljg (2002) investigated numer- 
ically also the dependence of the CF on fim and found it to be 
extremely weak. 

Defining the scaled velocity divergence, '& = Q^'^-d = 
— 57"" '' V ■ V, C04 found the following fit for the ratio of the power 
spectra: 

7^(fc) = [1 + (7.071fc)'']~" , (37) 



where 

a = -0.06574 + 0.29195cr8 for 0.3 < ag < 1 . 



(38) 



C04 argued that the ratio of the power spectra practically does not 
depend on the background cosmological model. This ratio is unity 
in the linear regime (k ^ 1) but decreases in the nonlinear regime, 
because the velocity grows slower than it would be expected from 
the linear approximation. 

In this paper we use Equations l l35t and i37l as the formulas 
respectively for the CF and the ratio of the power spectra. For ag we 
adopt the value obtained from a joint analysis of third-year WMAP 
and SDSS, ag = 0.77 (Spergel et al. 2007). In Equations i fTSt , 
i2U and i24\ , as the power spectrum we use a CDM spectrum. For 
zero baryon content, the shape parameter of the spectrum, F, equals 
simply to Qmh. Non-zero baryon content of the Universe modifies 
the shape parameter to (Sugiyama 1995) 

Fcfi = n^hexp [-fib (l + V2h/n,^)] . (39) 

Here we adopt Fcfr — 0.15, the value obtained both from first- 
year WMAP (Spergel et al. 2003) and a joint analysis of third-year 
WMAP and SDSS (Spergel et al. 2007). This value is in excel- 
lent agreement with the constraint on the shape of the power spec- 
trum of 2MASS galaxies, Fcff = 0.14 ± 0.02, obtained by Frith, 
Outram and Shanks (2005; assuming a flat ACDM cosmology, a 
primordial scale-invariant power spectrum and negligible neutrino 
mass). It is slightly higher than the corresponding result of Mailer 
et al. (2005), Feff = 0.12 ± 0.01, obtained using a measure of the 
three-dimensional power spectrum via an inversion of the 2MASS 
angular correlation function. 



4 GIfAVITY WINDOW OF 2MASS 

In this Section we derive the gravity window of the 2MASS survey. 
The 2MASS survey is dense, uniform and has an unprecedented 
sky coverage. Therefore, to a good accuracy it can be described 
by a spherical window. Since distances of 2MASS galaxies are un- 
known, the galaxies are weighted only by their fluxes and not by 
their distances (like, e.g., by the inverse of the selection function). 
This is the central assumption of the calculation below. 

The background light intensity due to uniform distribution of 
discrete sources is 



SdN, 



(40) 



where S = L/(47rr^) is the observed flux from the sources with 
intrinsic luminosity L and dN is the number of sources per stera- 
dian. In the case of uniformly distributed sources with the luminos- 
ity function (LF) $(!/), the contribution from a shell of thickness 
dr and radius r is dN — $(1/) dL r^dr. For a flux-, or magnitude- 
limited survey, only galaxies with L > 47rr^Smin are observed, 
where Smin is the limiting (minimal) flux. Hence, 

L 



e(r,L)- — =-$(L) dLr^dr, 



(41) 
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where 0{r, L) is the Heaviside step-function, QH(L — A-Kr^Sia\n)- 



Writing Lmi 

{L)No 



I = 



47r 



47rr Smin, this yields 
drWg(r). 



(42) 



(43) 



Here, iVo = di, and 

_ L$(L) di 

is the average luminosity of the population. The flux window of the 
survey is 

/^°° L^{L)dL 



(44) 



/„°°L<E.(L)dL ■ 

Wg gives the percentage of the total light from distance r which 
is included in the survey. Loosely speaking, it suppresses contribu- 
tions from distances larger than (L) / (47r5'min). Its detailed form 
is determined by the LF. 

The LF of 2MASS galaxies has been estimated by Bell et al. 
(2003), by matching a spectroscopic sample of Early Data Release 
SDSS galaxies with the 2MASS extended source catalog, to obtain 
redshifts for a subsample of 2MASS galaxies. Bell et al. fitted the 
2MASS LF by the Schechter function: 

$(L)dL = r (^)%xp(-i^)^, (45) 

where $* is the LF normalization, L* is the characteristic lumi- 
nosity at the 'knee' of the LF, where the form changes from expo- 
nential to power law, and a is the 'faint end slope'. For Ks-band, 
they found that a = —0.77 and the absolute magnitude A/*, cor- 
responding to the absolute luminosity L* , is M* = —23.29 + 
5 logiQ h. We adopt this form of the 2MASS LF here. 

The flux window of the 2MASS survey can be compared with 
the selection function of the survey, defined as 



(46) 



The selection function gives the probability that a randomly se- 
lected galaxy at distance r will be included in the survey. E06 call 
the selection function the 'number- weighted selection function', 
and the flux window the 'flux-weighted (or luminosity-weighted) 
selection function'. They note that "the number- weighted selec- 
tion function drops with the distance faster than the luminosity- 
weighted selection function. At large distances, we observe only 
the most luminous galaxies, so the amount of 'missing' luminos- 
ity from a volume of space is not as big as the number of 'miss- 
ing' galaxies". We fully agree. This implies in practice that when 
estimated from a galaxy survey, the flux dipole is a more robust 
quantity than the number dipole. 

If we want to exclude also the brightest sources, then 0(r, L) 
in Equation J41t becomes the product of two Heaviside functions, 
e(r,L) = eH(L-47rr25'„i„)-ejj(47rr-^Smax-L).Here,Smax 
is the upper limiting (maximal) flux. It is simple to check that the 
survey window then becomes 



Wg{r) = 



where L„ 



J^"-" L$(i) dL 
ffL<f{L)dL ' 



(47) 



E 47rr Smax. Thls windows suppresses also contri- 
butions from distances smaller than about -y/ {L)/(47rS'max), or, 
for the Schechter LF, just about L*/(47r5'max)- 



What would be the gravity window for the number dipole? 
The answer depends on the weighting scheme. In case of the num- 
ber dipole all galaxies are weighted equally - in a sense that they 
are not weighted by their fluxes or masses - but they may, or may 
not, be weighted proportionally to the inverse of the selection func- 
tion. If they are not, it is clear from the above analysis that then 
the gravity window is the selection function, 4>- (Now, instead of 
missing some percentage of the total light from distance r, we miss 
some percentage of all galaxies located there.) However, if they are 
weighted as l/(p{r), the gravity window is simply unity. This is 
so because the l/<j){r) weighting corrects for, on average, missing 
signal from large distances. The price to pay for this correction is 
huge variance of such an estimator of the LG gravity. This variance 
is called shot noise (see Subsection 15. U and has dominant contri- 
butions from large scales. Number-, rather than mass-, weighting 
of galaxies is another source of variance of this estimator, for sim- 
plicity also called shot noise. That shot noise comes predominantly 
from small scales. 

Using the IRAS 1.2 Jy survey, S92 measured the number 
dipole, weighting galaxies originally as l/(ji(r-). To mitigate shot 
noise and nonlinear effects from small scales and shot noise from 
large scales, S92 decided to modify these weights. They did this 
introducing the so-called standard IRAS window. 



WiRAS 




r < Ts , 

rs <r < R 

max ; 

R max < r . 



(48) 



This window is characterized by a small-scale smoothing, r^, and 
a sharp large-scale cutoff, 7?max. S99 adopted the values rs — 5 
h^^ Mpc and i?max = 150 Mpc, appropriate for the complete 
PSCz catalog (Saunders et al. 2000). Modified weights assigned to 
IRAS galaxies by S92 and S99 were 



weight (i) 



(49) 



It is clear that in this case, the gravity window of the IRAS number 
dipole is Wg. iras ~ Wiras- We will return to this point later. 

To specify completely the distribution for the misalignment 
angle (Eq.[27] or its small-angle approximation, Eq.|30t. we need 
the value of the velocity variance, Equation i2H . and of the cor- 
relation coefficient. Equation l l24t . In order to calculate the latter, 
we need to derive the Fourier form (Eq. [T6l l of the 2MASS win- 
dow, given above. In Equation l ll6t . we can use the fact that the 
spherical Bessel function, ji{x) — —{d/dx)jo{x), and integrate 
by parts. This yields 



(50) 



Wg{k) = / joikr)Wg{r)dr. 
Jo 

Let's cast Equation ( |47| ( to the form 



Wg{r) 



where 



J^--^{L)dL 



VE'(L) EE L$(L) . 

We can then write 

dWgir) dLn 



(51) 



(52) 



9I/max dr 

max ) 8nrSn 



dWgjr) dl/min 

dLmin dr 



(53) 
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Using Equations dSOt and ( I53t . and the fact that joix) — sinx/x, 
we obtain 



kJ^°°^{L)dLj, 



sin(fcr)'I'(Lmax)dr 
sin(fer)*(Lmi„)dr, (54) 



where ^{L) is defined by Equation \52\ . Note that the flux (or 
gravity) window does not appear in Equation MW for the velocity 
variance, and in Equation J24t for the correlation coefficient it ap- 
pears in such a way that its absolute normalization cancels out. In 
other words, the PDF for the misalignment angle is sensitive only 
to the shape of the gravity window. 

To relate the limiting fluxes to the limiting magnitudes, we re- 
mind that the observed minimal flux Smin is related to the apparent 
maximal magnitude i^max in the following way: 



SolO" 



(55) 



where 5*0 is the reference flux, which appears also in the relation 
between the absolute magnitude M* and absolute luminosity L* , 



M* = -2.51ogi 



47r(10pc)2So ■ 

In Equations (I55M56II we can therefore eliminate 5*0, obtaining 



(56) 



Smiu = 1.803 X 10" 



L* 



47r(l/!.-iMpc)2 



(57) 



Calculating the numerical coefficient in the above equation we have 
adopted M* = -23.29 + Slogm h (Bell et al. 2003). Follow- 
ing M03, for Kmax we have adopted the value 13.57. The reason 
for this choice of Kmax is twofold. First, our aim here is to im- 
prove the 2MASS window used by M03 properly accounting for 
nonlinear effects, which affect only the choice of optimal ifmin 
(the minimal magnitude). Second, M03 chose A'max = 13.57 be- 
cause "the extended source catalog is 97.5% complete within the 
SDSS early data release for extinction-corrected Rron magnitudes 
of Ks s; 13.57 mag" (Bell et al. 2003, Jarrett 2004). The 2MASS 
window is the 2MASS flux-weighted selection function under the 
assumption that the survey is complete within the flux limits. If we 
wanted to go deeper, we should account for increasing incomplete- 
ness as a function of distance. However, Figure 1 of M03, showing 
the convergence of the 2MASS dipole as a function of the limit- 
ing magnitude, suggests that contributions from al! galaxies (i.e. 
even from those not included in the survey) fainter than 13.57 mag 
are most likely negligible. Even for 2MASS galaxies brighter than 
13.57 mag (where the catalog is complete), "the faintest 300,000 
galaxies only change the dipole value by less than 5%" (M03). 
For Smax we have simply 



= S'n.inlO' 



0.4{13.57-i<'„i„) 



(58) 



If the brightest galaxies are not excluded, then either directly from 
Equation i44l . or from Equation l |54b , performing the limit Smax — * 
oo, we obtain 



8nSn 



kJ^<l>{L)dL 



sin(A:r-)*(Lmin)dr. 



(59) 



Figure [T] shows the 2MASS gravity windows for J^max = 
13.57 mag and different values of A'min. Dotted line corresponds 
to Equation ( |59t . i.e. to the case where the brightest galaxies are 
not excluded from the calculation of the flux dipole. Dashed and 




0.05 C.l 
k \h Mpc" 



Figure 1. Gravity windows for the 2MASS all-sky survey, for ii'max = 
13.57 mag and different values of ii'min- Dotted line corresponds to Equa- 
tion (59), i.e. to the case where the brightest galaxies are not excluded from 
the calculation of the flux dipole. Dashed and solid lines are plotted us- 
ing Equation (54j and describe the cases of excluding all 2MASS galaxies 



brighter than 



8 mae and K„ 



5 mag, respectively. Since 



all 2MASS windows have the same iCmax, they are similarly suppressed at 
large scales, i.e., at small k. For reference, we also plot the standard IRAS 
window (Eq. l6Qt . with Vs = 5 h ^ Mpc and Rmax = 150 ft- 1 Mpc 
(dot-long-dashed line). 



solid lines are plotted using Equation l l54t and describe respec- 
tively the cases of excluding all 2MASS galaxies brighter than 
Kmin ~ 8 mag (as done by M03), and -ftTmin ~ 5 mag (our 
choice, as justified below). Let's try to understand the influence 
of the limiting magnitudes on the shape of the gravity window. 
Since all 2MASS windows have the same i^max, they are simi- 
larly suppressed at large scales (small k). For a given A'min (corre- 
sponding to maximal limiting flux), all objects brighter than L* are 
excluded from distances smaller than r-min = L*/(47rS'max), 
and at distances r < rmin, all sources brighter than L*(r/rmin)^ 
are excluded. For the limiting magnitude Kmin = 5, rmin — 4.5 
ft" ^ Mpc, while for A'min = 8, rmin — 18.1 h'^Mpc. Conse- 
quently, the window for no exclusion of the brightest galaxies (dot- 
ted line) does not drop down at all for large k (small scales). The 
window for ifmin ~ 5 (solid line) does drop down but is fairly 
wide, while the window for A"min ~ 8 (dashed line) drops very 



rapidly. Since, as explained earlier. 
Equation l l48t we have 



Wg, IRAS = WiRAS, using 



iAs{k) = 



Sjijkrs) 
krs 



jo (kRn 



(60) 



For reference, we plot this standard IRAS window in Figure[T](dot- 
long-dashed line). The small scale smoothing of the IRAS window 
is rs — 5 Mpc, while for the 2MASS window with Kmin ~ 5, 
the effective smoothing scale rmin is about 4.5 /i~^Mpc. It is 
not therefore surprising that at small scales the IRAS window, ex- 
cept for its oscillatory behaviour, decreases fairly similarly to this 
2MASS window. 

Even neglecting shot noise, suppressing contributions to the 
flux dipole from small scales is necessary, since nonlinear effects 
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should be mitigated. For large k the coherence function of velocity 
with gravity (Eq. l35l l drops significantly below unity, decreasing 
the value of the cross-correlation coefficient (Eq. 1241 ). Then sup- 
pressing the gravity window for large k has almost no effect on the 
cross-term (which is the numerator of Eq. I24t . while it decreases 
the gravity variance, the square root of which appears in the denom- 
inator of this equation. This manipulation on the gravity window 
helps therefore to achieve the best possible correlation between the 
LG velocity and gravity. However, when one suppresses the gravity 
window for scales which are linear enough so that the CF is close 
to unity, one worsens the correlation again. This is so because even 
for linear fields (CF and the ratio of power spectra equal to unity) 
the correlation coefficient decreases for increasingly different win- 
dows of velocity and gravity. As a result, for some value of rmin, 
or Kmin, the correlation coefficient will have a maximum. 

We calculate the correlation coefficient. Equation i24\ (using 
the appropriate formulas for the CF and the ratio of the power 
spectra, and the velocity window given by Eqs. [T3] and I17t . for 
the 2MASS gravity window, for a range of values of the limiting 
magnitude -K'min. Results are shown in Figure |2] We see that p 
has a maximum (1 — p has a minimum) for Kmin — 4.5. Either 
not suppressing small scales at all (Kmin = — oo), or suppressing 
them excessively (Kmin greater than, say, 6) clearly decreases the 
correlation coefficient. This implies larger statistical errors of the 
estimated cosmological parameters when comparing the 2MASS 
dipole to the CMB dipole (see Sec.|6f; choosing the optimal value 
for -K'min is therefore very important. Instead of A'min = 4.5, as 
the optimal value we have adopted Kmin = 5. We have done this 
because the correlation coefficient changes in the range of Kmin 
from 4 to 5 hardly at all, while the number of excluded galaxies 
for Kmin < 5 would become very small, resulting in big Poisson 
noise. 

In our analysis so far we have not addressed the effect of 
shot noise. If shot noise is not negligible it also increases the opti- 
mal Kmin- This will be explained in Section |6] Shot noise for the 
2MASS flux dipole will be discussed in detail in Subsection 15. II 



5 OBSERVATIONAL ERRORS 

An estimate of the flux dipole from an all-sky galaxy survey is sub- 
ject to observational errors. One of them is shot-noise, due to di- 
lute sampling, by distant galaxies, of the underlying mass density 
field. Another one is the lack or deficit of galaxies in the Zone of 
Avoidance (at low Galactic latitudes). We will describe these errors 
correspondingly in Subsections l5.1l and l5.2l 

5.1 Shot noise 



The contribution to the LG gravity from a small volume element 
located at a distance r, estimated from the 2MASS survey is 




where a — T/(47rr^ denotes the sum over galaxies in- 

cluded in the survey and f is the unit vector towards the volume 
element. We model this quantity theoretically introducing a win- 
dow function, W: 



SM^—Y. Wiri)u,S^r^ ~ aW{r) ^ ) f . (62) 



0.07 I 




Figure 2. The con'elation coefficient, p (Eq. 1241 . for the 2MASS gravity 
window (Eq. l54L for a range of values of the limiting magnitude Kmin- On 
the ordinate, the value of 1 — p is plotted. The coefficient has a maximum 
(respectively, 1 — p has a minimum) for Kmin — 4.5. 

In the above we have accounted for the fact that the mass to light 
ratio for an individual galaxy, Ti, may not be equal to its average 
value, T, but may have some scatter. The quantity i^i = Ti/T; 
hence (vi) = 1. In Equation i62\ the summation is over all galaxies 
in the volume element, regardless whether or not they are in the 
2MASS survey 

The expectation value of the estimated gravity is 




= aWg{r){L)N{AV)r. (63) 

In the second step we have used the fact that the 2MASS gravity 
window we have constructed in Section|4l Wg, gives the percent- 
age of the total light from distance r which is included in the sur- 
vey. N{AV) is the number of all galaxies in the volume element 
(regardless whether or not included in the survey). The expectation 
value of the modelled gravity is 

(gA/) = aWir) ^(//,){i>)f = aWir){L)N{AV)r . (64) 

i 

In Equation i64\ we have assumed that the scatter in the mass- 
to-light ratio is independent of luminosity. However, the above re- 
sult for (gjv/) is also correct for fairly broad classes of luminosity- 
dependent scatter, e.g. for i/ = 1 + aJ-{L), where the luminosity- 
independent random variable a has zero mean. [In order for the 
variable f to be always positive, we have to impose an additional 
constraint |amin| < .?^max> where Qmin (negative since (a) = 0) 
is the minimum value of a and J-ma.x is the maximum value of the 
function JT.] 

Comparing Equation ( 1631 ) to ( |64t we see that if we adopt 
W = Wg, then the estimated gravity is an unbiased estimator of 
the modelled gravity. (From a different perspective, the modelled 
gravity is an unbiased estimator of what we really observe.) Our 
2MASS gravity window is constructed precisely in such a way to 
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assure this. However, the quantity g_E has scatter around gM- For 
flux (or number) dipoles estimated from flux-limited galaxy cat- 
alogs, the dilute sampling at large distances introduces significant 
scatter, called shot noise. The scatter in the values of g_B around gM 
due to the scatter in the mass-to-light ratio is not, properly speak- 
ing, shot noise (S92). However, for simplicity, we will call both 
these effects 'shot noise'. (S92 also follow this convention.) 

A lucid derivation of shot noise for the number dipole of IRAS 
galaxies can be found in Appendix A of S92. Shot noise for the 
2MASS dipole can be calculated in a similar way. However, except 
for the fact that the 2MASS dipole is a flux one, there is another 
important difference. While in the derivation of S92, galaxies are 
assigned weights essentially proportional to the inverse of the se- 
lection function, 2MASS galaxies are given equal weights (their 
distances are unknown). We will see below that this introduces a 
qualitative difference in the resulting formula for shot noise. 

To compute the variance of g_E we take the difference between 
g_B and gA/ for a full shell of thickness dr, we square it and calcu- 
late its expectation value. Finally, we sum up contributions to the 
total variance from all shells. The result is[f| 



2 

0"SN 



2M/g(n) + (l + Q) 



(65) 



Here, Q = {u'^) — 1^0 quantifies the amount of scatter in M/L. 
Were there no scatter, the value of Q would be zero. The function 
Fis 



F{t) = 



dL 



(66) 



The upper limit in the above integrals is either Lmax ~ 47rr^S'max 
or infinity, depending on whether we exclude the brightest objects 
or not (in the latter case, Smax = oo). 

Let us now investigate contributions to shot noise from small 
(r 0) and large (r oo) scales. For r ^ 0, consider first the 
case of no exclusion of the brightest objects. Then, both (j){r) and 
Wg(r) tend to unity. From Equation l l66t it is obvious that then 
-F(O) is a constant. Therefore, 



Ao-SN, nearby OC Q 



sf 



(67) 



nearby , i 



If Q is significantly greater than zero, then the RHS of the above 
proportionality blows up (since for 0, Si —> oo). This is 

shot noise from small scales, mentioned already in Section |4] It 
similarly plagues the number dipole (see Eq. 35 of S92, where there 
are similar r"** divergences). As already mentioned, to mitigate 
shot noise from small scales S92 introduced a window for the IRAS 
dipole. With inclusion of the IRAS window, contributions to shot 
noise from small scales in Equation (35) of S92 are proportional 
to WfjiAs{rt)/rf. For the standard IRAS window (Eq.|48ll, they 
scale as — » for 0, so shot noise from small scales is 

indeed strongly suppressed. As already stated, 2MASS galaxies are 
assigned equal weights. Still, shot noise from small scales can be 
mitigated. This is achieved by excluding from the calculation of the 
dipole contributions from the brightest objects, as described below. 
For finite S'max, the selection function is 



J-HL)dL ■ 



(68) 



Therefore, for r ^ and finite S'max, both Wg{r) and (f){r) tend to 
zero. However, although (f){r) tends to zero, it is straightforward to 
verify that the quantity W^{r) /4>{r) also tends to zero (at least for 
the Schechter form of the luminosity function). Finally, F{r) tends 
to a constant (though different from that for the case iS'max — oo). 
Hence, 



^(^In, nearby OC ^ ^ 



s 



(69) 



nearby , i 
Si<S,„ax 



The above sum is limited to objects with 5*; < S'max, what pre- 
vents it to blow up. Therefore, excluding the brightest objects is a 
good way to mitigate shot noise from small scales having at one's 
disposal angular data only. 

Contributions to shot noise from large scales do not depend 
on the choice whether we exclude the brightest objects, or not. For 
r —> 00, both 0(r) and Wg{r) tend to zero; it is straightforward to 
check that then also Wg{r)/(f>{r) tends to zero. The limit of F{r) 
for r — > oo is unity. Hence, using Equation i65\ we obtain 



distant 



St 



(70) 



We see that in the case of the flux dipole calculated with equal 
weights assigned to all galaxies, shot noise from large scales does 
not blow up; on the contrary, it decreases. As a result, in an anal- 
ysis of the 2MASS dipole one does not have to exclude any data 
from large distances. The analysis presented here assumed Kmin = 
13.57, but it is now clear that when calculating the 2MASS dipole 
one can include contributions from 2MASS galaxies fainter than 
this magnitude. (Although, as mentioned earlier. Fig. 1 of M03, 
showing the convergence of the dipole as a function of the limiting 
magnitude, suggests that their contribution will be negligible. See 
also Fig. 3 of Jarrett 2004.) 

Large-scale asymptotic behaviour of shot noise for the flux 
dipole (Eq. I70t with equal weighting is in contrast to the corre- 
sponding behaviour of shot noise for the number dipole calculated 
with galaxy weights proportional to the inverse of the selection 
function. In the latter case, contributions to shot noise from large 
scales diverge as <f)~^{ri) (see Eq. 35 of S92; at large distances 
(l>{rr) < 1)[!|To cure this problem, S92 introduced in their stan- 
dard IRAS window a sharp large-scale cutoff, i?max (see Eq. I48t . 
With inclusion of the IRAS window, contributions from large scales 
are proportional to Wfjijxs/4>'^{ri). IfWiRAS is truncated at some 
i?max, then there are no contributions to shot noise from scales be- 
yond i?max. Still, this does not imply that 'l/(^(r-)' weighting of 
distant galaxies is a good one. S92 were aware of this fact and con- 
structed the optimal window for the IRAS survey, i.e. such that it 
minimized variance of the estimator of the LG velocity (Eq. 45 
of S92). At large scales this window behaves asymptotically as 
J-j,{r)(t>{r), where J3(r) — JJ^ d'^r' ^(r') and ^ is the mass two- 
point correlation function. Therefore, the optimal weighting at large 
distances is proportional to J3 (r) (see Eq. 149 1 , so instead of in- 
creasing [as 4>~^{r) does] it decreases to zero, suppressing shot 
noise from large scales. (2MASS weighting is intermediate be- 
tween these two extremes.) At small scales, in the absence of the 



" In this paper we need a formula for shot noise only for illustrative pur- These divergences are due to the weighting scheme and not to the type of 

poses, therefore the derivation will be presented in follow-up work. the dipole. 
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scatter in the masses of galaxies, the window approaches unity. 
Therefore, the 1 /4>{r) weighting is then indeed the optimal one|j 
In the presence of scatter the window filters out small scales, as 
desired. Surprisingly, S92 resigned from using this window in the 
analysis of the LG acceleration and employed instead the standard 
IRAS window. The reason was that in the derivation of the opti- 
mal window they also attempted to account for nonlinear effects, 
but the coherence function they used was wrong (Chodorowski 
& Ciecieljg 2002). As a consequence, the resulting window fil- 
tered out small scales excessively. In the present paper, working 
with only angular data we have no choice: we have to assign equal 
weights to all galaxies. Though this is not the optimal weighting, 
this is still quite good: shot noise from large scales does not blow 
up. Moreover, excluding the brightest objects helps to mitigate shot 
noise from small scales. 

Let us recall: we denote the estimated gravity of the LG by g_B 
(Eq.|6It and its modelled gravity by g]\f (Eq.|62|). Although g_E is 
an unbiased estimator of gA/, it is still a biased estimator of the 
true gravity of the LG. Large depth of the 2MASS survey makes 
the estimated dipole to converge, but in order to mitigate shot noise 
and nonlinear effects from small scales we have to suppress con- 
tributions from small distances. (For angular data the only way to 
do this is to exclude the brightest galaxies, located preferentially 
nearby.) This reduction of the signal introduces bias in the estimate 
of the LG gravity. However, applying a Maximum Likelihood anal- 
ysis enables one to correct for this bias and to obtain an unbiased 
estimate of the parameter /3 — Q,^ /bL- This will be discussed in 
Section|6l 

How big is actual shot noise for the 2MASS survey? To an- 
swer this question, M03 performed bootstrap resampling on the 
2MASS galaxy catalog (100 times). They found that the standard 
deviation of the dipole direction was a fraction of a degree, and of 
the dipole magnitude a fraction of a percent. They concluded that 
'the systematic uncertainties are much larger than the shot noise'. 
Shot noise is certainly less an issue for the 2MASS dipole than 
for the IRAS PSCz dipole (S99) and for the 2MRS dipole (E06). It 
is smaller for the 2MASS dipole partly due to much bigger num- 
ber of galaxies in this survey compared to IRAS PSCz and 2MRS: 
there are about 13, 000 galaxies in the PSCz catalog and 23, 000 
galaxies in 2MRS, while for the limiting magnitude Ks — 13.57, 
the 2MASS catalog contains about 740,000 galaxies (M03). The 
main reason, however, is non-weighting of galaxies when calculat- 
ing the 2MASS dipole. (E06 weighted 2MRS galaxies inversely to 
the 'flux- weighted selection function', or, in our terminology, the 
gravity window, Wg). Still, M03 analysed shot noise including all 
(so also the brightest) 2MASS galaxies. Therefore, it is somewhat 
surprising that they did not found a trace of shot noise from small 
scales. A forthcoming paper of some of us (Bilicki & Chodorowski, 
in preparation) will be devoted to the optimal measurement of the 
2MASS dipole. We are planning to reexamine carefully the issue of 
shot noise there. Specifically, we are going to repeat the bootstrap 
resampling analysis and to compare its results to our analytical for- 
mula for shot noise. Equation l l65b . 

At first sight, it may seem surprising that Equation l |65l ) can 
be used in the case of only angular data, since radial functions 



° In a related paper, Feldman, Kaiser & Peacock (1994) constructed the op- 
timal estimator for the density power spectrum inferred from redshift sur- 
veys. They derived a formula for the optimal weighting of galaxies (Eq. 
2.3.4 of Feldman et al. 1994). For small r the optimal weight behaves like 
</)~^ (r), while for large r it approaches asymptotically unity. 



(j>{r), Wg{r) and F(r) appear in it. However, these functions are 
uniquely determined by specifying -ft'max (corresponding to Smin), 
KmiYi (corresponding to Smax) and the luminosity function of the 
2MASS galaxies. As described before, this luminosity function has 
been estimated e.g. by Bell et al. (2003). The only data employed 
in Equation J65b are fluxes. Si. (One also needs an estimate of Q, 
quantifying the amount of scatter in M/L.) Since we do not have 
these data at our disposal yet, for the rest of this paper we will ac- 
cept the claim of M03 that shot noise for the 2MASS flux dipole is 
negligible. 



5.2 The mask 

The source of the biggest systematic error in 2MASS remains the 
lack or deficit of galaxies in the Zone of Avoidance (at low Galactic 
latitudes). M03 masked the region of the ZoA, and repopulated it 
with 'synthetic galaxies'. In one method they cloned the sky above 
and below the masked region. In another method, they filled 'the 
masked region with randomly chosen galaxies such that it has the 
same surface density as the unmasked area'. The first method gave 
a dipole pointing towards I — 263°, b = 40°. The second method 
resulted in a dipole pointing towards / = 266°, b = 47°. M03 
adopted the mean of these two measurements as the best-fit dipole. 
However, the error bars they attributed to the mask-filling uncer- 
tainty were somewhat underestimated. We will return to this point 
later. 

The misalignment can be fully represented as a two- 
dimensional vector lying on the celestial sphere. In the absence 
of shot noise, the total misalignment is a vectorial sum of the 
cosmologically-originated misalignment Gc, described in Sec- 
tion|2l and the misalignment due to mask, 0m '- 

= 6^ + 6^. (71) 

We have G,n = {Al, Ab), where I and b are respectively the Galac- 
tic longitude and latitude. Under the simplest assumption, the distri- 
bution function for dm is a bivariate Gaussian of two uncorrelated 
variables of the same variance: 

piAl, Ab) = (27r)- V-^ exp . (72) 



The distribution for the modulus 6m ~ V AP + Afe^ results im- 
mediately from Equation J72t . It is a Rayleigh distribution (cf. 
Eq.m, 

HM-^exp(^-|L^. (73) 

Let us now invert the above reasoning and apply it to the variable 
6c. Since the distribution for 9c is (approximately) Rayleigh, the 
distribution for dc is (approximately) a bivariate Gaussian. The 
variable 9 is therefore a sum of two independent bivariate Gaus- 
sians, which itself is a bivariate Gaussian (of uncorrelated vari- 
ables). Hence, the variable 9 is Rayleigh-distributed, with the pa- 
rameter 

6? ^9l+a'^. (74) 

Here, 9^ = {6l^)/2, and cr^ = (0/2. 

The parameter 9, is defined by Equation ( I3II > and determined 
by the LG velocity variance (Eq.l2U and the correlation coefficient 
(Eg. 1 24b . Let us find an estimate for the mask variance a^. We have 
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(75) 



where Xi = [U — lY + [hi — 6)^, and (T, 6) are the means for 
the sample. As stated above, M03 study the effects of two differ- 
ent methods of 'repopulating' the masked regions with galaxies, so 
N = 2. Then X2 — xi, hence a — V2xi; the factor y/2 mustn't 
be neglected. This yields (in degrees) 



a ~ 5.4° 

(as opposed to 3.4°, or 4.5°, finally adopted by M03). 



(76) 



In the above likelihood, the 'data' are the measured values of the 
LG gravity and velocity, gm and v^, respectively, as well as firn, i.e. 
cosine of the misalignment angle. The model parameters ag and p' 
depend solely on Sv depends neither on P nor on 

We have written down the expression for the likelihood only 
for illustrative purposes. Therefore, for simplicity we will restrict 
our analysis to the case of given ■ Then, to find a maximum of 
the likelihood we calculate its partial derivative with respect to /3 
and equate it to zero. This yields the following equation: 



3(1 



/2no2 , VSp'fJ-mgmVm 



_3tv. 



= 0. 



(82) 



6 LIKELIHOOD FOR /3 

We mentioned in Section [T] that a comparison between the CMB 
dipole and the 2MASS flux dipole (the latter given by Equation |7]( 
can serve as a method to measure the parameter (3 — /bh- 
Of course, it cannot be done by naive equating of the two dipoles: 
such an estimate would be biased. Here we outline a likelihood 
estimation of (5 (for details see C04). 

In a Bayesian approach, one ascribes a priori equal probabil- 
ities to values of unknown parameters, which allows us to express 
their likelihood function, given v and g of the LG, via the proba- 
bility distribution function for v and g: 



£(param.) = p(v, g | param.^ 



(77) 



As the parameters to be estimated here we adopt /3 and b^; p is 
given by Equation l[8j. Theoretical quantities in this distribution are 
CTg, CTv, and p. Since now we account for observational errors, the 
variance of a single spatial component of measured gravity, (jg, is 
a sum of the ID cosmological component, o^ ,., and errors, e^/3 
(e denoting 3D errors, including shot noise and the mask). Here, 
gravity is inferred from a galaxian, rather than mass, density field. 
Therefore, a^ ,. — h\{g'^)/'i, where (g^) = (g • g) is given by 



Equation dlSb . To sum up. 



(78) 



Errors in the measured velocity of the LG are negligible com- 
pared to those in the gravity. The relation between the physical 
velocity, Vph, and the scaled velocity used in this paper, v, is 
Vph — f^m^v, hence ID velocity variance is 



(79) 



3 3 ' 

where s% = {v'^) = (v • v) is given by Equation Finally, 
errors in the estimate of the LG gravity do not affect the cross- 
correlation between the LG gravity and velocity, but increase the 
gravity variance. This has the effect of lowering the value of the 
cross-correlation coefficient. Specifically, 



1 + 



-1/2 



(80) 



where p is given by Equation J24t . 

From Equation l[8]l, the logarithmic likelihood for /3 and 6^ 
takes the form: 

\nC(fi,bL) = -31n(27r) - 31n [crg6i,Sv(l - p'f^'^] - 31n/3 
1 fgl. Zvl, 2^/3p' 



2(1 



+ 



(81) 



The LG gravity, inferred from the 2MASS survey, is tightly coupled 
to its velocity: 1 — p' ^ 1 and 1 — /im ^ 1- (See Table[T] ^obs = 
16° corresponds to p,^ — 0.96). At first approximation we can 
therefore assume p' — pm = 1, hence 



bhSv 3m 



«2) 



1/2 



5m 



(83) 



Thus, the estimate of j3 is not just the ratio of the LG velocity to 
its gravity: it is modified by nonlinear effects (which affect (v^) 
through the function TZ), different observational windows (which 
affect differently (g^) and (f^)), and observational errors. If all 
these factors are properly accounted for, then the estimate of /3 is 
unbiased. 

An optimal estimator is such that is not only unbiased but 
also has minimal variance. Expanding the logarithmic likelihood 
(Eq. I81t around its maximum up to second order in /3 enables one 
to find an estimator of the variance of p. In the strong-coupling 
regime (1 — p' < 1, 1 — pm < 1), it is 



(84) 



If errors are constant, i.e. they do not depend on -f^min, then a 
minimum of the variance corresponds to a maximum of the cross- 
correlation coefficient p. (The dependence of (g^) on /fmin is very 
weak.) Including higher-order corrections to the above formula 
does not change this fact. In the present paper, errors are indeed 
constant: the error due to the mask obviously does not depend on 
ifmin and shot noise is assumed to be negligible. The window func- 
tion of the 2MASS survey we have constructed here maximizes p 
(see Fig. |2j. This is why we call this window, under the assump- 
tion of negligible shot noise, optimal. It exactly corresponds to the 
minimal expectation value of the misalignment angle. 

As mentioned earlier, in follow-up work we will estimate shot 
noise ourselves. If we find that it is in fact not negligible, then it 
will influence the optimal value of Kmin- Shot noise as a function 
of Kuih-i monotonically decreases (see Eq.|69t. The factor 1 — p^, 
starting from the value of Kmin which maximizes p, monotonically 
increases (see Fig.|2]l. The interplay between these two opposing 
effects in Equation l l84t shifts the optimal ifmin (corresponding to 
a minimum of the variance of the estimator of P) to a larger value, 
compared to the case of negligible shot noise. 
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Table 1. Parameters of the distribution for the misalignment angle between 
the 2MASS and CMB dipoles, for various forms of the 2MASS gravity 
window. First column shows the limiting magnitude of the excluded bright- 
est galaxies, -R'min ■ Second column shows the correlation coefficient of the 
LG velocity and gravity, p, calculated according to Equation (24). Third col- 
umn shows the characteristic value of the misahgnment angle (in radians), 
0*, calculated using Equation )31t . Fourth column shows the correspond- 
ing value of the misalignment angle (in radians), including observational 
errors due to the mask, 9'^, calculated using Equation 174 1 . Fifth column 
shows the expectation value of the misalignment angle (in degrees), {9). 
Sixth and seventh columns show, respectively, the quantiles d\o and 9go (in 
degrees), defining the confidence intervals of 10 and 90% (for details see 
text). The last column shows the observed values of the misalignment angle 
(in degrees), obtained with and without exclusion of the brightest galaxies. 









W [°] 


eio n 


n 


fobs [°] 


- 0.951 


0.123 


0.155 


11.1 


4.1 


19.0 


16.0 


8 0.901 


0.183 


0.206 


14.8 


5.4 


25.3 


5.2 


5 0.969 


0.097 


0.135 


9.7 


3.6 


16.6 





0.08 - 




7 RESULTING DISTRIBUTIONS FOR THE 
MISALIGNMENT 

Figure[3]shows the resulting PDFs for the misahgnment angle, for 
various forms of the 2MASS gravity window. Like previously, dot- 
ted line corresponds to the case where the brightest galaxies are not 
excluded from the calculation of the flux dipole. The vertical stripe 
shows the value of the misalignment between the CMB dipole and 
the 2MASS flux dipole as calculated by M03, including all galaxies 
brighter than Ks — 13.57. The 'observed' value (16°) is greater 
than the expectation value for the angle (11.1°), but smaller than 
6lgo = 19.0° (see Table[B. 

To decrease the misalignment, in the second step M03 ex- 
cluded from the analysis all galaxies brighter than /("min = 8 
mag. A dashed line is plotted for the gravity window corresponding 
to this case. Consistently with Figure [2] the expectation value for 
the angle does not decrease; on the contrary, it increases to 14.8° 
(Table [T). Consequently, one would then expect the misalignment 
rather to increase. However, M03 noticed a substantial decrease of 
the misalignment, to about 5.2°. This value is smaller than the cor- 
responding 9io = 5.4°. Therefore, there is less than 10% chance 
that the decrease might have been accidental. Rather, an error in the 
analysis is more likely. Indeed, E06 repeated the procedure of MOB 
for the 2MRS data and essentially did not observe the decrease of 
the misalignment. 2MRS survey misses faint galaxies (fainter than 
Ks — 11.25), but does not miss bright galaxies. Therefore, if the 
effect was real, one should observe it also when using the 2MRS 
data. 

Solid line in Figure |3] is plotted using the window excluding 
2MASS galaxies brighter than Kmin = 5 mag. As described in 
Section 151 we expect this window to be close to optimal. Indeed, 
the resulting distribution is the narrowest among the three plotted; 
the expectation value of the misalignment drops to 9.7° and ^90 — 
16.6° (TablefTJ. Therefore, with (almost) 90% confidence we can 
expect the angle to decrease when performing such a preselection 
on 2MASS galaxies. Of course, this assumes constant mass-to-light 
ratio for all remaining (i.e., included) galaxies. 

The window with i^min ~ 8 mag is not optimal because it 
excessively mitigates nonlinear effects. This window excludes too 



Figure 3. Probability distribution functions for the misalignment angle be- 
tween the 2MASS and CMB dipoles, for various forms of the 2MASS 
gravity window. Dotted line is for the case where the brightest galaxies are 
not excluded from the calculation of the flux dipole. Dashed line describes 
the case of excluding all 2MASS galaxies brighter than ii'min = 8 mag; 
solid line describes the case of excluding all 2MASS galaxies brighter than 
^min = 5 mag. Thick solid vertical hne shows the measured value of 
the misalignment between the CMB dipole and the 2MASS flux dipole as 
calculated by M03, including all galaxies brighter than Ks = 13.57 (16°). 

many galaxies: while the window with K^^in ~ 5 mag excludes all 
Lt and brighter galaxies closer to the LG than about 5 Mpc, 
the window with A'min = 8 mag does the same for the distance 
of about 18 h^^ Mpc (see Sec. |4j. The signal from scales 5-18 

Mpc is sufficiently 'linear' to increase (if included) the corre- 
lation between the LG velocity and its measured gravity. 

M03 found 375 2MASS galaxies brighter than Amin = 

— 3/2 

8 mag. Based on this number and the relation N^^^ci oc Smax , 
where A'cxci is the number of excluded galaxies, we predict about 
six 2MASS galaxies to be brighter than Kmin = 5 mag. To reduce 
the misalignment calculated using their sample, E06 excluded five 
the brightest galaxies in 2MRS. They noticed a significant decrease 
of the misalignment, from about 21° to 14°. The five most lumi- 
nous galaxies in 2MRS are also the five most luminous galaxies in 
2MASS. Therefore, exclusion of these galaxies should work also 
for denser and deeper 2MASS survey. 



8 SUMMARY 

• An ultimate goal of comparing the CMB dipole to the 2MASS 
dipole is an estimation of the cosmological parameter /3 = 

• To obtain an unbiased estimate of (3, a good and standard 
method is Maximum Likelihood. 

• An important ingredient of this Likelihood analysis is the ob- 
servational window through which the 2MASS flux dipole is mea- 
sured, called here the gravity window of 2MASS. This window 
should be properly modelled. 

• By definition, the optimal window minimizes variance in the 
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estimate of optimizing the 2MASS window is therefore impor- 
tant. 

• In this paper, we have modelled the 2MASS gravity window 
and optimized it under the assumption of negligible shot noise. This 
optimization has been achieved by excluding contributions to the 
dipole from the brightest galaxies (Sec.©. Such an exclusion mit- 
igates nonlinear effects from small scales, which decorrelate the 
LG velocity from the estimated LG gravity. We have found that 
the optimal value of the minimal limiting magnitude, i^min (corre- 
sponding to maximal limiting flux), is about 5. 

• We have also demonstrated how to optimize the window in 
presence of shot noise. We have shown that the optimal value of 
K-mi-a will increase. 

• The misalignment angle is a sensitive measure of the correla- 
tion between the two dipoles: the higher the correlation, the smaller 
the expectation value of the angle (Eqs. |3H432| (. We have shown 
that a minimum of the misalignment corresponds to minimal vari- 
ance of the estimator of (3. A minimum of the misalignment is thus 
a sign of the optimal gravity window. 

• We have modelled analytically the probability distribution 
function for the misalignment angle (Sec. (2] App. [A^. We have 
shown that the misalignment estimated by M03 is consistent with 
the assumed underlying model (though it is greater than the expec- 
tation value). We have predicted that the misalignment is likely to 
decrease if 2MASS galaxies brighter than ifmin = 5 mag are ex- 
cluded from the calculation of the flux dipole. This prediction has 
been indirectly confirmed by the results of E06. 

• In a future work, we plan to perform the optimal measurement 
of the value of /3 by comparing the CMB dipole to the 2MASS 
dipole. We will thus have to fully specify the optimal window 
in presence of shot noise (though M03 claim that shot noise of 
2MASS survey is negligible). An estimate of shot noise can be ob- 
tained using methods described in Subsection 15.11 However, the 
misalignment angle can be used as an alternative way of optimiz- 
ing the window. As a function of A'min, the measured value of the 
misalignment will - with some scatter - initially decrease, reach 
a minimum and then increase (see Fig. |2]l. It is now clear that the 
value of ifmiii for which the measured misalignment has a mini- 
mum will be close to that optimizing the measurement of /3. 
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APPENDIX A: BEYOND THE SMALL-ANGLE LIMIT 

Here we check the accuracy of the small-angle approximation. 
Equation ( |30t , of the distribution function for the misalignment an- 
gle, given in general by Equation First, integrating by parts 
one can show that for v ^ 1, 



(Al) 



This yields for the error function a higher-order expansion (than 
Eq.[2l}: 



erf (s) ~ 1 - _ - -— s » 1 . (A2) 

\^ V^s 2 Vtts^ ) 

Using this expansion in Equation \T1\ for /i > 0, we obtain 

p(elv)=smeh+-^jexpf-^j, (A3) 

where 6, is given by Equation | |3U and Q < tt/2. For 6, 0, 
this distribution simplifies to the Rayleigh distribution (Eq.|30t, as 
expected. 

For ^ < 0, one can show in a similar way that 



erf(g^) 



-1 



f 1 



q > l.(A4) 



Using the latter expansion in Equation i27\ yields 

p(6llv) = for 7r/2 < e TT . (A5) 

Density distribution JA3b - iA5t has analytical cumulative distribu- 
tion function: 



F{e) 



1 — cos exp 
1, 



n/2 <e i^n. 



(A6) 



This allows for a straightforward estimation of the quantiles. More- 
over, using the fact that p = dF/dO and integrating by parts, the 
expectation value of the angle can be readily calculated: 



v) = yie.erf {e:') 



6** - o eU 



(A7) 



Thus, for Ol ^ 1, distribution JA3t -l lA5t has (almost) identical 
mean to the Rayleigh distribution (Ea.l32t. 

In practice, the total misalignment angle is a convolution of 
the cosmologically-originated misalignment, 6c, and the misalign- 
ment due to mask, 9m.. This convolution is not as simple as when 
both 6m. and 6c are bivariate Gaussians. However, non-Gaussianity 
of Oc + 6 m is smaller than of the variable 6'c with 0'^ = ^Jffl+cj^ 
(because 6m is Gaussian). We see in Table [T] that 9', is at most 
0.2. For el < 0.2, both distribution ([ASll-llASj and the Rayleigh 
distribution approximate the exact one (Eq. |27t very well. Specif- 
ically, since then < 1, we can expand distribution l |A3t -( |A5t , 
obtaining 



26^ 



where 



M{9) ^i + e? -le^ + +0 ("'^^ 

b bo/ 



(A8) 



(A9) 



It is straightforward to check that approximate distribution dASt is 
properly normalized. A simple calculation yields 



(AlO) 



Hence, for 9^ = 0.2, the Rayleigh distribution approximates the 
second moment of distribution JA8t (which, in turn, is then an ex- 
cellent approximation of the exact distribution) to 1.3% accuracy. 
Similarly simple calculations can be performed for other even mo- 
ments. 

Summing up, the distribution function for the misalignment 
angle between the CMB and 2MASS dipoles can be very well ap- 
proximated by its small-angle limit, Rayleigh distribution. It may 
be worth noting for other applications that, while for 6'^ > 0.2 the 
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exact distribution starts to deviate from the Rayleigli form, it is still 
well approximated by distribution ( IA3b - fA5l >. up to 9', ~ 0.5 (ex- 
cept for the very tail). In particular, the values of the quantiles ^lo 
and ^90 and of the mean angle remain within 2% from the exact 
values. 
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